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Abstract 

Evolution in finite populations is often modelled using the classical Moran pro¬ 
cess. Over the last ten years this methodology has been extended to structured 
populations using evolutionary graph theory. An important question in any such 
population, is whether a rare mutant has a higher or lower chance of fixating (the 
fixation probability) than the Moran probability, i.e. that from the original Moran 
model, which represents an unstructured population. As evolutionary graph theory 
has developed, different ways of considering the interactions between individuals 
through a graph and an associated matrix of weights have been considered, as have 
a number of important dynamics. In this paper we revisit the original paper on evo¬ 
lutionary graph theory in light of these extensions to consider these developments 
in an integrated way. In particular we find general criteria for when an evolutionary 
graph with general weights satisfies the Moran probability for the set of six common 
evolutionary dynamics. 


1 Introduction 


When modelling population evolution we are concerned with the spread of heritable 
characteristics in successive generations. The type of model that is used depends upon 
whether the population size is assumed to be finite or infinite. The majority of classical 
evolutionary models (see for example 0 , 0 ) use infinite populations, although finite 
population models are also well estiablished, the most important models being those 
in ii- These models are stochastic, and are solved using classical Markov chain 
methodology 0,B0. See also 0, 0 for an extension to evolutionary games in finite 
populations. 

The populations in the models described above, however, were “well-mixed”, i.e. 
every individual was equally likely to encounter every other individual. Real populations 
of course contain structural elements, such as geographical location or social relationship, 
which mean that some pairs individuals are more likely to interact than others. In such 
circumstances we need to be able to identify distinct individuals (or at least distinct 
classes of individuals), and considering finite populations is perhaps more natural than 
infinite ones (although finite structures each containing an infinite number of individuals, 
so called “island models”, were considered in [10). In 111 the modelling ideas of 0 


were extended to consider such structured populations based upon graphs, known as 
evolutionary graph theory. This has proved very successful, spawning a large number of 
papers (for example [3 El, 111 [isl . [la . 13 . [3 [10 ). For informative reviews see mmi. 

In an evolving population, we need to consider the mechanism of how the population 
changes, called the dynamics. Informally, the dynamics specify the way in which herita¬ 
ble characteristics are passed on from one generation to the next. For infinite populations 
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the classical replicator equation [22| is often used (although there are a number of alter¬ 
natives) , and in the stochastic model of Q there is a natural replacement dynamics built 
in. For structured populations this issue is actually considerably more complex, and the 
order of births and deaths, and where selection acts, is of vital importance [ 11 , 123 . We 


shall consider a set of dynamics that are commonly used in evolutionary graph theory 
models. The relationship between dynamics and structure is of key interest because the 
spread of heritable characteristics is directly dependent upon it. Whilst having essen¬ 
tially no effect on populations with no structure, this relationship potentially yields very 
different results on graphs. 

Under some circumstances it is, however, possible for the dynamics and structure 
to interact in such a way that the spread of heritable characteristics behaves just as if 
the population was homogeneous. This was a central theme of the classic paper ll|. 


where two important results, the circulation theorem and the isothermal theorem, were 
developed that addressed this question (see also for related work). In this paper we 
generalise the work of ll| to obtain a complete classification of when the combination 
of a population structure and dynamics can be regarded as equivalent to a homgeneous 
population in a precisely defined way, for the six most common evolutionary dynamics 
and graphs with general weights. 


2 The Model 


We shall first describe the population model of 0 , which generalises the model of iby 
incorporating a replacement structure. The notation used in this paper is summarised 
in Table [T] 

The population has a constant size N G h, N > 2, consisting of individuals 
Every individual is either of type A or B. 

This implies that there are 2^ different states of the population given by the combi¬ 
nation of type A and B individuals. We represent each state by a set S such that n £ S 
if an individual In is of type A. We can easily revert to using the number of type A 
individuals, l^l, if the population is homogeneous. The states 0 and Af = {1,2,. ..,iV} 
have only type B and A individuals respectively. 

Individuals have a constant fitness that may depend upon their type. 

The fitness of individuals in state S is thus given by the vector F(5') = (F'„(5'))„^^ ^ 
where 


,N 


Fn{S) = 


1 n ^ S, 

r £ (0, 00 ) n £ S, 


is the fitness of Here the fitness r of a type A individual is given relative to the fitness 
of a type B individual assumed to be 1. 

During a stochastic replacement event (that happens in an instant) an exact copy of an 
individual li replaces an individual Ij . 

The replacement events may be restricted in the sense that not all individuals can 
replace one another. To enforce such restrictions, [llj | imposed a replacement structure 
using a weighted directed graph given by the tuple {D, w) where D = (V, E) is a directed 
graph, with sets V of vertices and E of directed edges, and w is a map that assigns a 
weight to each edge such that w : V x V ^ [0, 00 ) : {i,j) >->■ Wij. Each vertex n £V 
represents therefore V = (1, 2,..., iV} so \V\ = N. We assume that {i, j) £ E ii and 
only if Wij > 0, which indicates that Ii can replace Ij. Note that we allow wu > 0 and 
therefore Ii can replace itself. All the information contained within the weighted digraph 
{D,w) is conveniently summarised by the N x N weighted adjacency matrix W = (wij) 
and therefore we will refer to {D,w) using W, which we call the replacement matrix. 

The replacement events are stochastic which means that there is a probability Xij = 
rij(F(5'),W) associated with (a copy of) A replacing Ij. There are several potential 
evolutionary dynamics on graphs that govern how the probability is determined. There 
three main types of dynamics that are summarised below, see also 2l|. We use the 
convention that L 


is chosen for birth and Ij is chosen for death. 

1. Birth-Death (BD): A is chosen first then Ij. We have that i £ 
probability bi and then (i,j) £ Ei is chosen with probability d. 




V is chosen with 
where Ei are all 
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Summary of Notation 


Symbol 

Definition 

Description 

N 

GZ+\{0,1} 

Population size. 

A,B 


The two types of individuals in population. 

In 


Individual n. 

S 

= {n: In of type A} 

State of the population. 

M 

= {1,2,...,7V} 

State in which all /„ of type A. 

r 

G (0, oo) 

Fitness of a type A individual. 

Fn{S) 

G {l,r} 

Fitness of In in state S. 

D 

= {V,E) 

Replacement digraph with vertices V where |y| = 
and directed edges E. 

Wij 

G [0, oo) 

Edge weight such that Wij > 0 if and only if 
(z, j) G E. 

W 

= (Wij) 

Replacement matrix: N x N weighted adjacency ma¬ 
trix of tuple (D,w). 

T+ 


Out temperature: Sum of all outgoing edge weights 
of vertex n € V. 

T- 

= Ei=l W^n 

In temperature: Sum of all incoming edge weights of 
vertex n € V. 

h 

G [0, 1] 

Probability li chosen for birth. 

dij 

G [0, 1] 

Probability a copy of R replaces Ij given R was cho¬ 
sen for birth, i.e. replacement by death. 

dj 

G [0, 1] 

Probability R chosen for death. 

bij 

G [0, 1] 

Probability a copy of R replaces R given R is chosen 
for death, i.e. replacement by birth. 


G [0, 1] 

Probability a copy of R replaces R. 

Pss' 

G [0, 1] 

State transition probability. 

S 

= {Pss') 

State transition matrix. 

^=(=,W,r 


Stochastic process with state transition matrix S 
such that * dynamics are used on graph W and type 
A individuals have fitness r. 

Ps 

G [0, 1] 

Fixation probability of type A individual given initial 
state S. 

W 


Set of all strongly connected replacement matrices. 

Wc 

{W : T+ = T- Vn} 

Replacement matrices that are circulations. 

Wi 

(W : 7;+ = T- Vz,j} 

Replacement matrices that are isothermal. 

Wh 

{W : r+ = 1 Vn} 

Right stochastic replacement matrices. 

Wl 

{W : T- = 1 Vn} 

Left stochastic replacement matrices. 

Cn 


Replacement matrices whose digraphs are cycles of 
length N. 

fn 

(Wij) ^ {Wij/J2n'^in) 

Map from W to Wr. 

h 

(Wij) ^ {Wij/Y,nWnj) 

Map from W to Wr. 

f 

(Wij) ^ {Wij/J2n,kWnk) 

Map from IF to IF. 

M, 

Replacement matrices for which £1* is p-equivalent to 
a Moran process when * dynamics are used. 


Table 1: Notation used in this paper. 


edges starting in vertex i. dij is used to signify that there is ‘replacement by death’. 
Finally, = bAj. 

2. Death-Birth (DB): Ij is chosen first then We have that j & V \s chosen with 
probability dj and then (i,j) S Ej is chosen with probability bij, where Ej are all 
edges ending in vertex j. bij is used to signify that there is ‘replacement by birth’. 
Finally, = dib^j. 

3. Link (L): li and Ij are chosen simultaneously. In this case {i,j) G E is simply 
chosen with probability tij. 

For each type of these dynamics, the natural selection can, through the fitness parameter. 
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Process 

F(Ii replaces Ij) 

Order chosen 

P(Chosen first) 

P(Chosen second) 

BDB [11] 

Zij — hi dij 

li then Ij 

— WS) 

w,, 

‘ E^"(^) 

> 

BDD [2^ 

Xij — hidij 

li then Ij 

n 

bi = — 

n 

^ _ w^j/Fj{S) 

'' Y^W^n/FniS) 

DBD [IJ 

II 

Ij then li 


n 

b- - 



DBB [IJ 

II 

Ij then li 

n 

dj = T7 

TV 

n 

_ WijFi{S) 

'' Y^njF^iS) 

LB [IL] 

^ WijFi{S) 

Simultaneous 

N/A 

n 

N/A 


LD [2^ 

n,k 

^ _ w^j/Fj{S) 

Simultaneous 

N/A 

N/A 

y^,Wnk/Fk{S) 

n,k 


Table 2: List of the dynamics used in this paper. Note that L will be used in place of 
LB and LD where appropriate. 


influence either the choice at birth (resulting in adding “B”) or at death (adding “D”). 
It yields 6 kinds of evolutionary dynamics on graphs summarized in Table [5J These 
dynamics have been extensively studied, in particular, see 123 f or a detailed comparison 
of them. Of these, the BDB and LB dynamics were used in |ll| . 


2.1 The fixation probability 

The fixation probability, pg = pg(*, W,r), is the probability that the population with 
initial state S is absorbed in J\f where * is the dynamics being used. 

Given that the replacement events are random, the transitions between the states of 
the population are described by a stochastic process, which we denote £. The properties 
of £ can be investigated once the state transition probabilities of moving from state S to 
S', Pss' = are calculated using the replacement probabilities as follows: 


Pss' = 


ies 

E ’ 

ijes 


(F(S), W) if S" = S' \ {j} for some j G S, 
(F(S), W) if S" = S U {j} for some j ^ S, 
Cy(F(S),W) iiS' = S. 


The transition probabilities, Pss', satisfy the Markov property because they only depend 
upon the state S, that is, the probability of transitioning from the present state to another 
state is independent of any past and future state of the population. The stochastic process 
£*.w.r with state transition matrix S = S(=i=,W,r) = iPsS')s,S'c{i,2,...,N} is therefore 
a Markov chain. The Markov chain S*,w,r is part of the class of evolutionary Markov 


chains described in [27 1. 


The absorbing states of are 0,A/", which means that if the population is in 

either one of these states then it remains there indefinitely. This property of can 

be used to measure the success of a type A individual by calculating the probability that 
it fixates, that is, everyone in the population is of type A. The fixation probability is 
then given by solving 


Ps = 


E 

S'C{1,2... 


n A 

PsS'Ps' 


( 1 ) 


,N} 
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with boundary conditions Og = 0 and pff = 1 . 

As demonstrated in [2^ , LB and LD dynamics may differ in time scale but they yield 
the same fixation probabilities when fitness is constant (which is our case). Thus, for our 
purposes the dynamics are the same and we will thus consider them together and denote 
them by L. 

2.2 The Moran Process 

The Moran process Q , a stochastic birth-death process on finite fixed homogenous pop¬ 
ulation, can be reconstructed as fBDB.WH.r- for a constant replacement matrix 

Wh = (1/iV),,,. (2) 

For any r € (0, oo) and any S C {1,..., N}, the fixation probability for this process, or 
Moran probability, is given by 


r 1 — r 

Ps = I 1 - r-^ 

i|5|/fV 


iir^l, 
if r = 1 . 


We are interested in characterizing graphs (and evolutionary dynamics) that yield the 
same fixation probabilities as the homogeneous matrix Wh given in ([2]). We note that for 
this matrix all of the transition probabilities Vij take the same value independent of i,j 
or the dynamics, and consequently the fixation probability under each of the dynamics 
is the same. 


2.3 Classes of Graphs/ Matrices 

The set of all admissible replacement matrices is defined as follows 


W = {W : for every i,j, there is n such that (W"')ij > 0 }. 


This definition means that W is strongly connected as for any pair of vertices i and j, 
there is a path (of length n) going from j to j. Unless specihed otherwise, we will consider 
admissible replacement matrices only. 

As in lll|, for any W (admissible or not) we dehne the in temperature of In, T~, and 


the out temperature of T+, by 


N N 

i=i i=i 

W is called a circulation if , for all n G V and it is called isothermal if 

T/" = T~, for all i,j G V. W is called right stochastic if T+ = 1, for all n gV and it 
is called left stochastic if T~ = 1, for all n G V. The sets of all circulations, isothermal 
matrices, right stochastic matrices, and left stochastic matrices, respectively are denoted 
by Wc,Wi,Wr, and Wl respectively. 

The set Cn denotes the sets of matrices representing cycles of length N, more specif¬ 
ically, for (wij) G Cn we have wu = 1/2 for i = 1,2,.. .N, = ■ ■ ■ = = ■ ■ ■ = 

■u^iAT-iiN = WiNii = 1/2 for some permutation U, * 2 ) ■ • ■, of the sequence 1 , 2 ,..., A^, 
and Wij = 0 otherwise. 

We also define the maps Jr'.W ^ Wr, /l : W ^ Wr, and f : W ^ W respectively, 
by 




( Wjj \ 

\ Srt / 


fh {{Wij)) 


f Wjj \ 

VE„ Wnj) 


and f'{{w^j)) = 


. Sn.fc 


Note that Jr preserves right stochastic matrices and Jr preserves left stochastic matrices. 
Moreover, fR^W) — /l(W) for all W G Wr. Also, since /' simply involves multiplying 
W by the constant 1/ f. Wnk, it implies that W G Wc ^ //W) G Wc- 
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When the dynamics *, matrices Wi and W 2 , and fitness r are given, we say that an 
evolutionary Markov chain is p-equivalent to f*^w 2 ,r if for every S C {1,..., N}, 

Pg{*,Wx,r) = pg{*,W 2 ,r), in which case we write Wi W 2 . 

We are specifically interested in finding matrices equivalent to the Moran process. 
For a dynamics *, we define 

M* = {W : W Wh for all r > 0}. 


3 Results 

The map fn preserves the equivalence classes of BDB and BDD dynamics, /l preserves 
the equivalence classes of DBB and DBD dynamics and /' preserves the equivalence 
classes for link dynamics. Specifically, as one can see from the proofs in the Appendix, 
for any W and any r > 0 

W ^bdb,. /ij(W), (3) 

W ^BDD.r /ij(W), 

W -DBB.r /l(W), 

W ^dbd,. /l(W), 
w /'(W). 

We thus obtain the following results, which completely specify the graphs which are 
equivalent to the homogeneous matrix Wh for each of our evolutionary dynamics. 

Proposition 1 (Link). = Wc- More precisely, the following statements are equiv¬ 
alent: 

(a) W is a circulation. 

(b) For all r > 0, W '^L,r Wh. 

(c) There is r > 0 such that W Wh. 

We note that Wc = f'~^{Wc) = {W; /'(W) G Wc} and thus, similarly to Proposi¬ 
tion [5] below, Proposition [1] can be written as = f'~^(Wc). 

Proposition 2 (BDB and DBD). Mbdb = and Mdbd = /^^(l^c)- More 

precisely, the following statements are equivalent: 

(a) fnfW) is a circulation. 

(h) For all r > 0, W ^bdb.t 

(c) There is r > 0 such that W ^bdb.t Wh 

The equivalent conditions for DBD are similar to the above for BDB but f^ is replaced 

by h- 

Proposition 3 (BDD and DBB). Mbdb = //{^({WhIuCat) and Mdbb = ^({Wh}U 

Cn) . More precisely, the following statements are equivalent: 

(a) fniW) = Wh or /«(W) G Cn- 

(b) For all r > 0, W ^BDD,r ^Ah. 

The equivalent conditions for DBB are similar to the above for BDD but fB is replaced 

by h- 

In particular, Mbdd C Mbdb and Mdbb C Mdbd- The sets M, are illustrated in Figure 

m 

Note that unlike in Propositions [1] and [H Proposition [3] does not contain “any r 
implies all r”. In fact, when r = 1, there is no selection and thus the dynamics BDB and 
BDD are the same (and also the dynamics DBB and DBD are the same). Consequently, 
by Proposition [21 

W ^bdd.i Wh ^ /fl(W) G Wc o W G Mbdb, 

W ^dbb.i Wh o /l(W) g Wc o W g Mdbd- 
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w 


Ws 


Wb 

Wi 

W4 

W7 

W2 

Ws 


W3 



Wi = Wi nu Cjv) 

= M//n/^^({W h} u Cjv) 

W2 = W/\/^l({WH}UC'iv) 

= m \ u c^) 

Wi = Wc\ Wi 

W4 = {fn^Wc) \ Wc) n /«H{Wh} U Cn) 
W5 = (fu^Wc) \ Wc) \ /fi'({WH} U Cm) 
W "6 = (fE^Wc) \ Wc) n /^^{Wh} u c^) 
Wr = (fr^Wc) \ Wc) \ /i ^({Wh} U Cm) 
W"8 = M/ \ ULi Wi 


We 

Wi 

W4 

W7 

W2 

W5 


W3 



Ml = Wc 



Wi 

Wi 

W7 

W2 

W5 


Ws 



Mobd 


Ws 






We 

Wi 

W4 


W7 

W2 

W5 



W3 



Mbdb 

Ws 






We 

Wi 

W4 



W7 

W2 

Ws 




W3 









Mbbb 



Mbbb 


Wg 


We 

Wi 

W4 

W7 

W2 

Ws 


W3 



Mbbb n Mdbb 


Figure 1: The diagram on top shows eight partitions Wi of W (labelled using their index 
i = 1, 2,..., 8). A combination of these partitions make up the sets M*. Below this, a 
separate diagram for each of the standard dynamics is given showing the partitions that 
make up M*. The bottom right diagram shows the partition where £ is p-equivalent to 
a Moran process regardless of the standard dynamics on the graph being used, that is. 
Ml n Mbdb n Mbdd H Mdbd H Mdbb = Mbdd H Mdbb- 


3.1 Our results in the context of known results 


For the LB dynamics, Proposition [T] was stated and proved in [llj as the Circulation 


theorem. For the LD dynamics, Proposition [T] follows from the Circulation theorem and 
the result of [ 2 ^ that the fixation probabilities for LB and LD are the same. 

As shown in Appendix (AT] BDB is the same as the LB dynamics for right stochastic 
matrices (in particular, for BDB dynamics. Proposition [5] can be seen as the Isothermal 
theorem from [HI). Proposition [5] thus follows from Proposition [I] thanks to ([21). The 
natural symmetries between fn and /l and BDB and DBD dynamics allow us to extend 
the Isothermal theorem to DBD dynamics as well (see also [28|). 

Overall, Propositions [1] and 5] and the occurrence of Wc within them are consis¬ 
tent with the claim made in [11 1 that the circulation criterion completely classifies all 
replacement matrices where is p-equivalent to a Moran process. 

Our most important new result is Proposition |21 It shows that the BDD and DBB 
dynamics require very strict conditions to yield the Moran process. Either the population 
structure is homogeneous, or it is a directed cycle. This latter structure is an interesting 
theoretical example, but is unlikely to apply to real populations, meaning that the ho¬ 
mogeneous population is practically the only way to get the Moran process for a realistic 
population. 
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3.2 The importance of self-loops in BDD and DBB dynamics 

Proposition [3] by definition requires that wu > 0 Vi = 1,2,... ,7V. Without such self¬ 
loops, fBDD.w.r, ^DBB,w,r cannot ever be p-equivalent to the Moran process. The ability 
of an individual to replace itself therefore plays an important role in the replacement 
structure of the population and cannot be discounted. For BD dynamics, when increasing 
the diagonal weights of W, the fixation probability decreases for BDB and increases for 
BDD. For DB dynamics, the increase in fixation probability DBB is greater than that 
for DBD. For LB dynamics, the fixation probability remains the same. 

With BDD and DBD evolutionary dynamics on graphs one may encounter the fol¬ 
lowing problems if there are no self-loops [l^, page 245]. For DBB dynamics, a type A 
individual with almost infinite fitness still has a fixation probability bounded away from 
1 because even type A individuals can be randomly picked for death and replaced by 
type B individuals. With self-loops, however, a type A individual will almost always be 
replaced by itself (or another type A individual) and therefore has a fixation probability 
approaching 1. Similarly, for BDD dynamics, a type A individual with almost zero fitness 
does not have near probability 0 of fixating as type A individuals can be randomly picked 
for birth and replace type B individuals. With self-loops, such an individual will almost 
always pick itself (or another type A) to replace and therefore its fixation probability is 
near 0. Thus the inclusion of self-loops removes some problematic features of the BDD 
and DBB dynamics, and makes them more attractive dynamics to use in models. 


4 Discussion 


In this paper we have considered an evolutionary graph theory model of a population 
involving general weights and a variety of evolutionary dynamics based upon the work 
of 111, which was a development of the classical population model of Q. In such pop¬ 
ulations, the population size is fixed at all times and at successive discrete time points 
one replacement event occurs. Like the aforementioned papers we consider two types 
of individuals, where fitness depends upon type but no other factors (i.e. there are no 
game-theoretic interactions). In particular the single most important property of such a 
process is the fixation probability, the probability that a randomly placed mutant indi¬ 
vidual of the second type will eventually completely replace the population of the first 
type. 

This fixation probability depends upon the fitnesses of the two types of individuals, 
but can also be heavily influenced by the population structure as given by the weights, 
and by the evolutionary dynamics used. These effects are commonly observed, although 
in some circumstances evolution proceeds as if as on a well-mixed population as from 
the original work of Q, dependent only upon the fitnesses of the two types, and some 
important results in this regard were already given in llj . The aim of this paper was to 
provide a generalised set of conditions for when this would be the case. 

By dehning what is meant by fixation-equivalence to the Moran process, we provided a 
general result which, independent of the specific dynamics used, helps identify graphs that 
do not affect the fixation probability. With respect to each of the standard dynamics, 
we then classified sets of evolutionary graphs that have the same fixation probability 
as the Moran process (or well mixed population). These sets include graphs that are 
circulations and therefore generalises the work of [HI- 

An important new result shows that the set of weights for which we obtain fixation 
equivalence to the Moran process for the BDD and DBB dynamics is very restricted, and 
so that for most populations with any structure this equivalence will not hold for these 
dynamics. We note also that the inclusion of non-zero self weights wa eliminates some 
problematic features of these two dynamics (i.e. that individuals with 0 htness could 
fixate or those with infinite htness could be eliminated) and so improves the applicability 
of these dynamics. 

Presenting evolutionary dynamics on graphs in the way that we have allows one to 
incorporate a variety of dynamics in their analysis, both of standard type and other 
dehnitions. This will improve our understanding of dynamics on graphs in general. We 
note that the list of dynamics in Table [2] is not exhaustive. For example, [s^ used 




imitation dynamics, which is a class of DBB dynamics with an additional requirement 

Wii > 0 Vd 

In general the inclusion of non-zero self weights, in contrast to many earlier evolution¬ 
ary graph theory works, allows for a greater flexibility of modelling. We note that this 
is consistent with the original work of Q, which allowed self-replacement as an integral 
part of the process. For well-mixed populations it does not matter much whether this 
possibility is included or not (at least for sufficiently large populations with intermediate 
htness values), and it is likely that it has often been excluded for reasons of convenience 
because of this without the ramifications being fully considered in many later works. It 
is thus important to consider whether to include such self weights when modelling spatial 
structure using evolutionary graph theory. 
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Appendix 
A Proofs 


A.l BDB is the same as LB for right stochastic matrices 


For BDB dynamics we have = bidij. By definition bidij = 1, we can therefore 
write this as = b^d^j jY,n,k bndn,k ■ Substituting b^ = Fi j J2m=i gives 

d F F 


d F- 

LLij ± I 


'^n,k {dnkFji yX]m=l 


Er, 


> djik^n 


If W is right stochastic, i.e. 1 ^ for BDB dynamics 

we have that dij = Wij giving = WijFi Ij2n,k WnkFn which is the 

LB dynamics as required. We also have that DBD is the same as LD for left stochastic 
matrices. The explanation follows the same procedure as above. 
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A.2 Lemma [T] (Forward Bias) 

The key Lemma [T] stated below is used in the proofs of all propositions and it relies 
heavily on the notion of forward bias of state S which is then given by the ratio of the 
probabilities of a forward transition to a backward transition from S. A forward and 
backward transition from S occurs when the number of type A individuals increase and 
decrease by one respectively, which happen with probability 


-Ps = X! ^s,Su{n} and Pg Ps,s\{n}- 

n^S tiGS 

Lemma 1 (Constant Forward Bias). Let £ he an evolutionary process on states S C 
{1, 2,..., iV} with transition probabilities Ps,S' that satisfy 

• Ps,S’ > 0 only if S and S' differ in at most one element 

• for every S 0, |1,..., iV|, there are S'^ and S~ such that IS’+I = IS’I + 1 and 

\S-\ = \S\-l andPs,s+>O.Ps,S->Q■ 

Then, the following are equivalent 

a) There is a constant c > 0 such that for all S' C {1, 2,, N} 


Ps = \ 1 - c-^ 
[\S\/N 


ifc^l, 
if c = I 


h) £ has constant forward bias, that is, there is a constant d such that for all S C 
{1,2,...,1V} 

P+ /p- = d. 

Moreover, if either (a) or (b) hold, then c = d. 

Note that a similar result is given in where the forward bias is explicitly 

defined as 

aGSb^S / aeSb^S 

which is what one gets when using Link dynamics, or BDB dynamics if W € Wr. Note 
that in Lemma [T] the forward bias is defined independent of the dynamics and therefore 
applies to all dynamics that satisfy the assumptions. 

Proof, “(a) (b)”: Take any S C (1, 2,..., iV}. It is known that 


Ps = 


Ps,S'Ps' - Ps,SPs + X! (^ASUlnlPsuln}) + X! {Ps,S\{n}Ps\{n} 

S' n^S tiGS 

D-\- 


and using Ps^s = 1 — Pg — Pg gives 

0 = X! (^Ps,Su{n} ~ Ps')') +^2 (^^S,S\{n} (^Ps\{n} “ PS 

n£S 


For c 1, equation (g]) simplifies to 
1 — — 1 + 


(4) 


0 = 


1 — c 


-N 


PtP 




1 — c 


.-N 


PI 


p+ /p- — 1 _ 

MS - „-isi- 


:-|S| _ C-l^l + l 1-c 


C-ISI-1 _c-|S| c-1-1 

For c = 1, equation (g]) simplifies to 

0 = i\S\ + 1 - \S\)P+ + (|,S| - 1 - \S\)Ps => P+ /Pg = 1. 
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“(b) (a)”: The state transition matrix S = {Ps,s') can be scaled to give S' = {P's s') 

such that T5 s = 0 and P's^s' — Ps,S’ /(I — Ps,s) = Ps,S’/{Ps + ) where S' is a non¬ 

absorbing state. The fixation probability Ps will be the same whether S' or S is used. 
This is because equation m can be rearranged as follows 

Ps = X! ^ss'Ps' ^ Ps= PssPs + X! Pss’Ps' 

S' s'-.s'^s 

Ps{^-Pss)= X! Pss'Ps'^ Ps= X! p+ Tp- Ps>- 
s'-.s'^s s'-.s'^s -s ‘S' 

Let {So,Si,... ,S7v} be a partition of the states S such that S S Si if |S| = i. The 
probability Pij (S) of transitioning from state S € Si to lumped state Sj with respect to 
S' is 


P^AS) 


f0 j ^i±l, 

j = i — 1, for i = 1, 2,..., TV — 1. 
[d/{d+l) j = i + l 


This can be easily verified, for example, take j = i — 1 then 


P.,_i(S)= ^ P's,s> 

S'eSi-i 


E 


Ps,S' 


Ps + Ps 


Ps 1 

-f Pg 1 + d 


(5) 


since the forward bias is equal to d. Equation satisfies the necessary and sufficient 
condition for the Markov chain with state transition matrix S' to be lumpable with 
respect to the partition {So, Si,... ,Sn} (Theorem 6.3.2 page 124, [HI). Let S = (P,j ) 
be the state transition matrix for this lumped Markov chain then the probability Pij of 
transitioning from lumped states Si to Sj is given by 

Pi,j = Pi,j{S)- 


The state transition matrix S describes a random walk with absorbing barriers and 
therefore the probability pf of type A individuals fixating when the population starts in 
lumped state Si is calculated using the methods in Q to give 


Pi 


i-l j 

1+En 




Pk,k—1 

Pk,k+1 


N-1 J 


i+En 

i=l k=l 


Pk,k—1 

Pk,k+1 


In this case. 


1 -d"* 
1 -d-^ 
i/N 


d=l 


since Pk._k-i/Pk,k+i = 1/r for A: = 1, 2,..., TV — 1. By definition, pg = p^ where i = |S| 
as required. □ 


A.3 Proposition [T] (Link) 

The following statements are equivalent: 

(a) W is a circulation. 

(b) For all r > 0, W -L.r Wh- 

(c) There is r > 0 such that W ^L,r Wh. 

(d) For all r > 0 and for all S C {1, 2,..., TV}, the forward bias of fn.w.r is r, i.e. 

Ps/Ps =r. 
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(e) There is r > 0 such that for all a € {1, 2,..., A^}, the forward bias of the one 
element set S = {a} is r, i.e. 


^{a},{a,b} 

b^a 


Pa,0 


= r. 


Proof. For LB dynamics the forward bias is given by 


p+ 

Ps 


WabPa 

aes b^S L/ '^'^kFn 


WbaPb 

aeS b^S L/ ^’^kPn 

n.k 


aGSb^S 
aesb^S 


For LD dynamics the forward bias is given by 


p+ 

Ps 


EE 

aesb^S 


Wab/Fb 

^ ^ '^nk/Fk 
n,k 


EE 

aes b^S 


Wba/Fg 

^ ^ '^nk /Fk 


n,k 


aeSbt^S 

EE'^'>“ 

oes b^s 


“(a) ^ (d)”: W is a circulation i.e. T+ = T„ for all n e {1,..., N} and thus 
aes b^S aes ^ n kGS ' aGS ^ kGS ' 

P,'Y,Wab = Y^ (Pa~-J^Wka) = f W„a - Wka) => 

aGSb^S aGS ^ kGS ^ aGS ^ n kGS ^ 

'P,'Y,Wab = P^'Y^Wba- 

aGSb^S aGSb^S 

Note that ^ because W is admissible and represents a strongly con¬ 

nected graph. Thus, the forward bias for both LB and LD is equal to r. 

“(d)=>(e)” is trivial as (d) is much stronger than (e). 

“(e)=>(a)” Let a and r is fixed. By above calculations of the forward bias, we have 


N N 

Wab= y] Wba ^ -Waa + = -Wga + '^Wia=^ 

&^S'—{a} {a} i—1 i—1 

therefore W is a circulation. 

“(d)=^>(b)” follows from Lemma[TJ 
“(b)^(c)” is trivial. 

“(c)^(e)” follows from Lemma [T] 


N N 

y] Wai = y] W^a 

i=l i=l 


□ 


A.4 Proposition (BDB and DBD) 

More precisely, the following statements are equivalent: 

(a) /fl:(W) is a circulation. 

(b) For all r > 0, W ~BDB.r Wr. 

(c) There is r > 0 such that W ^sDB.r Wr 

(d) For all r > 0 and for all S C {1, 2,..., iV}, the forward bias of fBOB.w.r is r, i.e. 

Ps/Ps =r- 
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(e) There is r > 0 such that for all a G {1,2,..., iV}, the forward bias of feDB.w.r of 
the one element set S = {a} is r, i.e. 

P{a},{a,b} 

b^a 

- 

Pa3 


Proof. Let U = {uij) = /fl;(W) = [wij/ then for BDB dynamics the forward 

bias oi SBDB.MV.r is given by 


p+ 

Ls_ 

Ps 


EE 

aeSb^S 


Fg Wab 

^ ^ Fn ^ ^ ’^an 


EE 

aeSbfS 


Fb Wba 

^ ^ Fn ^ ^ 


E E 

ggS b(S 

EE“b“ 

b^S aeS 


and therefore the forward bias of ^BDB.w.r is the same as forward bias of £^BDB,u,r- 
Similarly, with almost identical working as above, when V = /l(W), the forward 
bias of foBD.w.r is the same as forward bias of foBO.v.r and is given by 


p+ 

Ps 


EE 

aes bfS 


1/Ft Wab 

^ ^ 1/ Fn ^ ^ 'l^nb 


EE 

a^Sb^S 


n n 

l/i^g Wba 

^ ^ l/-^n ^ ^ '^na 
n n 




a&Sb^S 


^EE^-' 

aes b^S 


and the proof of the Proposition for DBD closely follows the one for BDB given below 
with U and fn appropriately replaced by V and Jl- 

“(a)=> (d)”: If U = /i?(W) G ITc, i.e. if U is doubly stochastic, then the forward 
bias (for S ^ 0,^/") is equal to 


El 

Ps 


Cl^S 


'E ( 


fees 


aeS ^ n kes 



" E E 

a£Sk^S 


i-s”! - EE^fc“ 

aeskes 


= r 


“(d)=>(e)” is trivial as (d) is stronger than (e). 

“(e)^(a)” Let a and r is fixed. By above calculations of the forward bias, we have 


EE''“''=EE''''“' 

aeSbfS aeSb^S 


Consider the states S' = {a} in which there is only one individual of type A then 


N N N 

^ ^ '^ab — ^ ^ y^ba ^ '^aa F ^ ^ y^ai — Uga F ^ ^ Uig 1 = ^ ^ Uig 

b^S b^S i=l i=l i=l 

is true for all a = 1,2,..., and therefore U is doubly stochastic and thus /fl:(W) is a 
circulation. 

“(d)=J>(b)” follows from Lemma[TJ 
“(b)^(c)” is trivial. 

“(c)=J'(e)” follows from Lemma [T] □ 
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A.5 Proposition [3] (BDD and DBB) 

The following statements are equivalent: 

(a) /fl(W) = Wh or /^(W) e Cn- 

(b) For all r > 0, W ~BDD.r Wr. 

Proof. The replacement probabilities (F(S'), W) for BDD dynamics can be rewritten as 
ry (F(S'),U) where U = (u^) = /fl:(W) = (wij/ by multiplying the numerator 

and denominator with Wm as follows 


r.,(F(5),W) = 

u^J/Fj{S) 
'Yhn '^in/Fn{S) 


W^3/Fj{S) 


_^ 1 Wjj/ {Fj{S)Y,n'^^n) 

N Y.n Win/Fn{S) N Win/{Fn{S) Yn 

= r.,(F(5),U) 


and therefore we have that W ^^BDD.r U, for all r > 0. The forward bias using U for 
state S is given by 


p+ 

Pc 


1 ^ab/Fb 

n 

1 Uba/Fa 

^SHS^'E^bn/Fn 


'^ab 


EEv ,F 

aeSb^S / .'^an/rn 


■EE 

a^Sb^S 


“^ba 


( 6 ) 


^ ^ '^bn/Fji 


Similarly, let V = (uy) = /l(W) = {wij/ Yn Then for DBB dynamics we have 

WijF^ _ WijFi/Yn'^nj _ VijFi 


'^n'^rtjFn 'l2nFn / Yn 

and therefore the forward bias when using 'V is given by 

1 VabF, 


Sn ^njFn 


p+ 


EE 

aeS b<^S 


Vab 




'•EEy- 

aeSbi^S / .’^nbi'n 


^ EE 

aesb^S 


EE' 


Vba 


1 '^ba^^b 

^ ^ VnaFn aeS b^S E 


The proof of the Proposition for DBB closely follows the one for BDD given below with 
U and fji appropriately replaced by V and /l- 

A.5.1 If U S C]\[, then U ~BDD,r 

If U S Cn then there are only two nonzero elements in each row. In particular, in row i 
of U we have that = 1/2 for some ki ^ i. In the numerator of equation (jS]) for 

a G S, b ^ S and ka ^ a we have that for all S 


'^ab 


'^ab 


0 if 6 ka, 

1/2 


Ugn / Fn{S) Uaa/Fa{S) + Uaka/Fk^{S) | l/2r+l/2 if^— ^a- 


Similarly, in the denominator of equation ([5]) for a G S , b ^ S and kb ^ b we have that 
for all S 


'^ba 


^ba 


0 if a 7^ kb, 

1/2 


Ubn/Fn{S) Ubb/Fb{S) + Ubkt, / Fki,(,S) l/2+l/2r if « — 
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This means that equation (jS]) for all S can be written as 

i/2r+-l/2/ rl/2+l/2r 

where x (y) is the number of nonzero Uab (uba) terms in the numerator (denominator). 
If we partition the vertices of the digraph of U into any two sets Vi , V2 then the number 
of edges e{i,j) and e{j,i) for i G Vi and j G V2 are by definition the same because it is 
a cycle. This means that for o G S' and b ^ S the number of nonzero Uab^Uba terms in 
the numerator and denominator respectively are the same hence x = y and rx/y = r as 
required. As per Lemma [H Sbdd,u,t- is p-equivalent to the Moran process. 


A.5.2 If U ^BDD.r- Wh for all r > 0, then U = Wh or U G Cn 

By Lemma [H the forward bias ([S]) is equal to r for all S C {1,..., N} giving 

Uab _ Uba ^ 

aeSh^S ggs b(f-S Ubn/Fn 


aes 2^ '> 
j'^s 


^ ^ ^ab 
b^S 


igS 


= E 


Uba 

aeS 


b^S ^ ^ Ubj 


1 


7E 

ies 


(7) 


'^bi 


Note that if r = 1, © holds for all U G Wc- From now, we will consider r 7^ 1 only. 
For clarity, the remainder of this section of the proof is broken down into the following 
six steps. 


Step 1: Derivation of general state dependent row-sum equation 

Let {7(a, S) = i-®- ^ — U{a, S) — J2j^sUaj- Equation © thus becomes 


1 — U (u, S) U (6, S) 

^ 1 - U{a, S) + Uia, S)/r " ^ 1 - Uib, S) + U{b, S)/r 

1 U{n,S) 

1 + U{a, S)(l/r - 1) = 1 + U{n, S)(l/r - 1) ’ 

Equation ([5]) can be written as a Taylor series as follows 


CXD N OC 

- 1)'= [U{a, 5)]'= = y] C/(n, 5) y](-l)yi/r - 1)^= [U{n, 5)]^= 

aGS k—0 n—1 k—0 

CXD N 00 


1 [U{a, S)f = YY(^- 1/0'= {U{n, 5)]'=+' 

aGSk—0 n—1k—0 


( 8 ) 


(9) 


For equation m to hold for all r the coefficients of (1 — 1/r)^ should be same, that is, 
for all k 


N 

Y{U{a,S)f = Y{U{n,St+\ 

a^S n—1 


( 10 ) 


Step 2: The diagonal of U consists of non-zero elements 

Consider the state S' = {a} then equation gives 

N 

uL^Y^nt^- (11) 

n—1 

If Uaa = 0 or 1 then m implies that all off-diagonal terms in column n are zero which 
is a contradiction with W (and thus also U = /fl(W)) being strongly connected, which 
means that 0 < Uaa < 1- 
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Step 3: The column of U contains nonzero elements, all equal to l/run 
Since 0 < Uaa < 1 , we can divide equation (ttH) by u^a giving 


We have that 



( 12 ) 


and therefore (fT^ implies that 0 < Una < Uaa- There must he n ^ a such that Una = Uaa 
as otherwise, by (HU), we would have Uaa = 1 - Let Ca = {i ■ Uia = Uaa }■ (HU) becomes 


1 = 


E 

ieCa 


E 




= \Ca\Ua 


E 

j^Ca 


,k+l 


( 13 ) 


As k ^ oo, (IT^ implies that Uaa = l/|Ca|- Thus, again by (|T^ . uja = 0 for all j ^ Ca- 
This means that in column n of U there should be = | C„ | with 2 < m„ < N nonzero 
elements, including that are all equal to l/m„. 


Step 4: m„ is the same for all n 

Considering state S = {i,j} and using Uaa = (fTPl) can be written as follows 

1 1 / 1 1 \ 

(u., + + {u,, + ^ j 

where a, / 3 ,7 are the number of rows where l/wi is adjacent to 0, 0 is adjacent to l/m-j, 
and l/wi is adjacent to l/m^ in columns i and j respectively. More precisely, a is the 
cardinality of the set = {n : Uni = llrru.Unj = 0}, fi is the cardinality of the set 
Kfj = {n : Uni = 0 ,Unj = 1/wj} and 7 is the cardinality of the set = {n : Uni = 
ljmi,Unj = l/wj}- 

Since Ci = U and C,- = Kij- U we have that ruj = a + 7 and m-j = (3 + 7. 
Since are disjoint, we have a + ,5 + 7 < N. Now, consider the different 

possibilities we can have on the left-hand side of equation m- 

Case 1: 

Uii = Ijuii, Uij = 0 in row i and Uji = llrrii, Ujj = l/rrij in row j. Thus 0,7 > 1 and 
therefore equation m gives 


-b 


rrii -b rrij 


i <2 *' 




+ 




m 


k-\-l 


rriimj 


fc+i 


a + ^ + 27 
(a + 7)^ ' V + 7)(/3 + 7) 


a 


(a+ 7 ) 


k-\-l 


iP + i) 


fc+1 


a -b /3 -b 27 \ 

^ 7« + 7)(/3 + 7)y 


(/3 -b 7)^" -b (a -b /3 -b 27)^" a{f 3 -b 7)^^+^ -b / 3 (a -b 7)'"+^ -b 7(0 -b /3 -b 27)^"+^ 


[(a-b7)(/3-b7)]'= 

(/I+ 7 )"+ (« + /?+ 27 )" = 


[(a-b7)(/3-b7)]'=+i 


k , t , « , o + / 3 (a + 7)^^+^ -b 7(0 + 13 + 27)'=+^ 


(/3 + 7)" + (a + /3 + 27)" = 


(a-b7)(^-b7) 

a (/5 + 7 )^ / 3 (a + 7 )* (07-b/ 37 -b 27^)(a-b/ 3 -b 27)^" 


a -b 7 


/3 + 7 


a /3 -b 07 -b ,57 -b 7^ 


liP + iY 

a -b 7 


/ 3 {a + 7)'= (72 - a/ 3 ){a +/3 + 27)'^ 


13 + 1 


al 3 -b 07 -b ,07 -b 7^ 


As fc —>■ 00, we get (/3 -b 7)^ (a -b 7)^ + {a + (3 + 27)^ since a-b/3-b27>/3-b7, a-b7 
hence we want 7^ = a /3 to get rid off {a + (3 + 27)^. This implies that ,5-b7 = a-b7=> 
a = /3 ^ a = [3 = 1 giving rrii = rrij. 
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Case 2: 

Uii = llrrii^ Uij = l/rrij in row i and Uji = 0, Ujj = ^/mj in row j. This case is 
symmetrical to Case 1 and therefore we get that a = /3 = 7 giving rrii = mj. 

Case 3: 

Uii = ^laii, Uij = l/rUj in row i and uji = Ijrrii, ujj = l/nij in row j. Thus 7 > 2 and 
therefore equation (HU gives 


rui + rrij 

miTrij 


a + /3 + 27 

(a + 7)(/3 + 7) 


a 




^k+l 


+ 7 


rrii + rrij 

miTrij 


fc+i 


+ 7)^+^ + / 3 (a + 7)^+1 + 7(a + /3 + 27)^+^ 
[(a + 7)(/3 + 7)]'=+i 


, a , ^ a{P + -if+^ + + +-i{a +l3 + 2-^f+^ 

2 (a + p + 27) = - 


2(„ + /J + 27)'-=“<'’ + ^' 


(a + 7)(/3 + 7) 

/ 3 (a + 7)'" ^ ^7 + /37 + 27^)(a + /3 + 27)'^ 


a + 7 /3 + 7 ' Q!/3 + 07 + /37 + 7^ 

(2a/3 + 07 + ^7)(a + ^ + 27)'= a(/3 + 7)'= /3(a + 7)'‘ 


aP + (37 + /37 + 7^ 


a + 7 


^ + 7 


As A: ^ 00, we get (a +/3 + 27)^ 7^ (^ + 7)^ + (a + 7)^ since a + /3 + 27 >/3 + 7, a + 7 
hence we want 2ap + a7 + /37 = 0 =7 a, /3 = 0 giving rrii = aij. 

Case 4- 

Uii = ^iTTii, Uij = 0 in row i and Uji = 0, Ujj = l/rrij in row j. Thus a, P > 1 and 
therefore equation (HI gives 


1 /to^ + l/mj 


a 


m. 


fe+i 


P ( mj + rrij \ 
V mirrij J 


1 1 

(a + 7)^ (/3 + 7)^ 

{P + 7 )* + (a + 7 )'' 

[(a + 7)(/3 + 7)]'= 

(/3 + 7)'= + (« + 7)fe: 

(/3 + 7)fc + (a + 7)'=: 


- Q; P / 7 + /3 + 27 

(a+ 7)'"+! (/3 + 7 )'‘+^ \(q; + 7)(/3 + 7), 

a (/5 + 7)^+^ + / 3 (a + 7)'"+! + 7(a + /3 + 27)'"+^ 

[(a + 7)(/3 + 7 )]'=+i 

«(/5 + 7)^+^ + P{a + 7)'=+! + 7(a + /3 + 27)'"+^ 

(a + 7)(/3 + 7) 

a (/3 + 7)^ _l_/ 3 {a + 7)^" 7 (a +/3 + 27)'"+^ 

q; + 7 /3 + 7 Qf /3 + 0:7 + /37 + 7^ 


fe+i 

=7 

=7 


As k ^ 00 , we get 0 {a + P + 27)^ since a,P > 1 hence we require that 7 = 0 to get 
an equality. 

Conclusion from all the cases above 

We see that ^ mj is potentially possible only in Case 4. However, U is strongly 
connected. If one connects i and j by a path i = iq, ii, *2, • ■ • *„ = j, then one has 
mi^. = mif.^-^ as ik and ik+i must fall into Case 1, Case 2 or Case 3 above. Thus 
mi = mj. This implies that every column of U has 2 < m < N nonzero elements, 
including u„„, that are all equal to 1/m. This is also true for every row of U because it 
is right stochastic by definition. 


Step 5: There exists state S such that Ca = Ca’ for all a,a' G S 

We can define the state TZx = {n ■ Uxn = Uxx] then, by definition, x G TZx and \Ti-x \ = m 
since there are m nonzero elements in row x of U. Consider the state S = TZx \ {y} for 
y G TZx \ {a:}- For this S (as well as any other state), we have that 


if n G S' then 
if n ^ S then 



m 
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We can therefore write equation (fTUl) in the form 


min(m,|S'|) 

2=1 



min(m,|5|) . . . fc+i 

2=0 ^ ' 


(15) 


where Xsii) is the number of U{n,S) terms equal to ijm for n G S and X'g{i) is the 
number of U{n,S) terms equal to i/m for n G M, which means that X'g(i) > Xs(i) for 
i 7^ 0. The ratio of the left-hand side and right-hand side of equation (ITSl) should always 
be equal to one. Therefore, as fc ^ oo, we require that 


As(*„,ax) = A's(^„ax) — 

TO 

where imax is the largest i such that Xs{i) > 0. 

We have that imax = to — 1 in equation m because [S'! = to — 1 so U{x,S) = 
(to — 1)/to. This means that for state S', as fc —>■ oo, we require that 

777 — 1 

Xs{m - 1) = X'sim - 1)-. 

m 

Since Xs{m — 1) is an integer, X'g{m — 1) has to be a multiple of to and the only possible 
value that satisfies this criteria is X'g(rn — 1) = m hence Xs{m — 1) = to — 1. 

Since Ag(TO—1) = to there exist to rows ji, J2, ■ ■ ■ ,jm such that U{jn, S) = (to—1)/to, 
that is, Uj^a = l/m\/a G S. This means that Ca = {ji, , Jm} Va G S hence Ca = Ca' 
for all a, a' G S. 


Step 6: m = 2 or m = N 

By contradiction, assume that 2 < m < N. We can consider another state S' = TZx \ {z} 
such that z G TZx \ {x,y}- We then have that imax = to — 1 in equation m because 
|S'| = TO — 1 so U{x,S') = (to — 1)/to. As before, this means that Ca = Ca' for all 
a, a' G S'. Since x G S,S' and TZx = S U S' we have that Ca = Ca' for all a, a' G TZx- For 
2 < m < N this implies that vertices i G TZx are disconnected from j G N \ TZx and we 
therefore have disconnected graphs, a contradiction. □ 
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